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Blocking is considered for an N-trunk group of exponential servers 
with Poisson-offered load whose rate parameter varies with time. The 
infinite trunk case is solved by means of a rapidly convergent series of 
Poisson-Charlier polynomials. This solution is used to obtain practical 
approximations of blocking probability, transition probabilities, and 
recovery function for general time-variable offered load in the finite trunk- 
group case. An integral equation is derived satisfied by the blocking 
probability in the general case. In the situation of constant offered load, 
two additional methods are derived for providing easily computable 
approximations; one based on the integral equation, the other based on an 
approximate inversion of the Laplace transform. To aid in the latter 
approximation, bounds on the roots of Poisson-Charlier polynomials are 
obtained; in particular, an approximation is obtained for the dominant 
root. The inversion of the integral equation is studied with the purpose of 
providing the basis for future investigations of errors of approximation. 
Curves are provided for a number of examples permitting comparison of 
exact and approximate solutio?is. 

I. INTRODUCTION 

The main purpose of this paper is to present a discussion of the 
time behavior of blocking in a fully available iV-trunk group for any 
initial state with exponential servers and with Poisson-offered load 
whose rate parameter a(t) itself may be considered to vary with time; 
that is, the probability of j calls arriving in a time interval (0, t) is as- 
sumed given by 

/ f ' \\ L a{u)du \ 
exp {- j o a(u)duj Ly ° ^ , j = 0,1,2, ■■ -. (1) 

The service rate is taken equal to 1, so that a(t) is measured in erlangs. 
The problem of blocking with time-variable offered load was con- 
sidered by Palm 1 for finite trunk groups and by Khintchine 2 for infinite 
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trunk groups. The impetus for this is the need felt for more accurate 
computation of blocking probabilities and correlation information 3 than 
can be obtained by quasi-stationary analyses, that is, by the use of 
equilibrium formulas in which the offered load parameter is replaced 
by its instantaneous value. Lack of statistical equilibrium renders 
this approach inaccurate. The time-variable aspect of the input stream 
should be carefully distinguished from other statistical descriptions 
such as peakedness, 4 since the effects on the system are separately 
identifiable. It has been reported, for example, that offered load and 
peakedness determinations from carried usage, peg count attempts 
offered to a group, and overflows are misled by the time variability 
of the offered load. 

Palm had proposed an interesting method of accounting for the time 
variability of the offered load, i.e., his "slow variations" model. In 
this model, it was assumed that the actual ordering in time of a(t), 
that is, the functional dependence of a(t) on t could be ignored if a(t) 
varied slowly. He replaced a(t) by a random variable with an incom- 
plete gamma function distribution. 1 Thus, traffic functions such as 
blocking may be obtained from their equilibrium values by averaging 
over the appropriate gamma distribution. This model, however, re- 
quires further elaboration in view of the investigations of Iversen, 3 
who showed that the Palm approach does not correctly model the 
empirical data collected in the extensive Holbaek measurements of 
Danish telephone traffic. Iversen found that the correct time variation 
of the traffic could not be ignored. 

The trunk provisioning procedure whereby one uses the average 
offered load over a busy hour to achieve a required grade of service 
results, in some cases that were considered, in only a small under- 
estimate of the required number of trunks as calculated by the methods 
of this paper. Since the standard method is convenient, this may be 
viewed as substantiation of the approach. 

Essential for the methods of this paper is a Volterra integral equa- 
tion derived in Section II satisfied by the blocking probability, 
Pff(t, N), experienced by a load of a(t) erlangs offered to an iV-trunk 
group. Exact analytical solution of this equation is not useful, but 
numerical methods may be advantageously used. An important fea- 
ture of the equation, nonetheless, is that it permits studying errors of 
approximation and, in one instance (Appendix D), was directly used 
in the construction of an approximate operator for studying the 
transient response in the case of constant a(t). Appendix A presents an 
explicit representation of PN(t, N) for general a(t) by means of an 
infinite Neumann expansion. Inequalities for Pxit, N) and truncation 
error estimates for the Neumann series are also given. 
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The infinite trunk group case forms the basis of the approximations 
developed in Section IV that are applicable to the general case of time 
variable a(t). Although this case was solved by Khintchine, 2 Section III 
presents new representations in terms of rapidly convergent Poisson- 
Charlier expansions. Truncation error estimates are obtained, and the 
rate at which the state probabilities approach the Poisson form is 
assessed. To aid in the use of Poisson-Charlier polynomials, Appendix 
B provides a short discussion of their properties, especially providing 
convenient means of expanding a function into a Poisson-Charlier 
series. Since the state probabilities of the infinite trunk group system 
are often close to Poisson, this form of representation is very useful. 
The Poisson-Charlier expansion expresses the deviation of a function 
from the Poisson form. Further, in Section IV, the Poisson-Charlier 
polynomials are used to express the transition probabilities in explicit, 
closed form. 

The approximations of Section IV are applicable to time variable 
a(t), and are developed from the infinite trunk group solution by 
renormalization appropriate to the finite trunk group. To facilitate the 
use of the approximations, closed expressions are obtained for the 
relevant infinite trunk group solutions. This approximation procedure 
gives rise to the useful notion of a "modified offered load." One of the 
approximations obtained was, in fact, already obtained by Palm. 1 
This approximation is particularly interesting because it uses the 
Erlang loss function, B (N, a), for which rapid methods of computation 
are available. 8 - 6 A special case of the approximations for transition 
probabilities is that for the recovery function, which is important in 
the discussion of correlation properties 7 and, hence, in the determina- 
tion of variances of traffic parameter estimators. 

The constant offered load case is studied in Appendix C. Although 
the solution for the state probabilities is known, 8 the integral equation 
formulation appears to be new. Certain advantages are obtainable from 
this formulation. The errors of approximations to the state probabilities 
satisfy the same integral equation but with a different inhomogeneous 
term; thus, the more general integral equation is studied, leading to 
methods for investigating the quality of approximation. For this 
purpose, a natural Banach space (uniform norm over [0, <»]) is 
introduced, in which the integral operator is bounded and has a 
bounded inverse. Of course, the known Laplace transform of the 
transition probabilities is immediately obtained as a corollary. The 
integral equation is also used in Appendix D as a tool for the con- 
struction of an approximate solution (the scaling method) correspond- 
ing to an arbitrary initial state. Appendix C also presents several 
bounds on the required roots of Poisson-Charlier polynomials by 
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obtaining general bounds on the largest and smallest members of sets 
of positive numbers subject to isoperimetric constraints; the results 
are then specialized to the Poisson-Charlier polynomials. One of the 
bounds for the dominant root is explicit and easily calculable. Its 
accuracy appears to be good (see Fig. 8). 

Appendix D provides two approximations for the case of constant 
offered load. The scaling approximation, whose genesis was suggested 
by S. Horing, is constructed by obtaining an approximate time in- 
variant of the transition blocking probability from the initially empty 
state. The subsequent generalization of this approximation to arbi- 
trary initial states is then obtained by means of the integral equation. 
The Laplace transform approximation is constructed by an adaptation 
of the Widder formula 9 for the inversion of the Laplace transform. It 
requires the determination of the dominant root; but, depending on 
the needs, it may be made arbitrarily accurate. It has the interesting 
property that, under certain conditions, it provides bounds for the 
exact solution. 

A discussion of some results and graphical illustrations is given in 
Section V. In testing the quality of the modified offered load approxi- 
mations, high change rates of a(t) were chosen, in fact, much higher 
than would occur in practice. The errors of approximation increase with 
increasing rate of change of a(t), hence, the examples chosen indicate 
much higher errors than one would expect to encounter in the practical 
application of these methods. 

II. INTEGRAL EQUATION FOR BLOCKING 

It is the object of this section to establish Theorems 1 and 2, which 
provide integral representations of the binomial moments (13) of the 
probability distribution (2) of the number of busy trunks, and corollary 
1 of Theorem 2, in which an integral equation is given for the prob- 
ability that all trunks are busy at a given time. 

Let £ = £(t, N) be the number of trunks busy at time t in an iV-trunk 
group, and P, = Pj(t, N) the corresponding probability, 

Plkit, N) = n = P&, N). (2) 

The probability generating function g(t, f, N) is given by 

g = m Ht.N) = £ P .( t} JV)^. (3) 

3=0 

At the point t + dt, one has 

g(t + dt, f, N) = g + dg, |(* + dt, N) = £ + d£ (4) 
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and, hence, 

g+dg = Et{tSEK*i\lQ} J (5) 

in which E[£ d * | £] is the conditional expectation of f d * given £, and E$ 
is the expectation over the probability distribution of £. The boundary 
at j = N necessitates a further analysis of (5). One has 

g + dg = (1 - PN)Et[t*Ett d t\t < iV]} 

+ Pjr*itt*SQ*»|$- O. (6) 
The probability distribution of eZ£ is 

P[d£ - - 1 1 $ = i, £ j < #] = jdt, 
PLdl = 0|£ = j, ^ j < Nl = 1 - (a + j)<ft, 
P[<** = 1 1 € = i, ^ i < JT| = adt, (7) 

p[d£ = o | £ = JV] = 1 - Ndt, 
P[d£ = - 1|£ = AT] = iVdi; 
hence, 

^Cf d « I « = ;, £ j < AT] = 1 + (of - a - j + jr 1 )*, 

i?Cf*«|{ = iv] - i - jv(i - r l )dt. { } 

From (6) and (8), one has 

ff = (1 - PN)E£!*(at - a - * + &- 1 ) | { < AT] 



- PtftfiCff*" 1 ^ - 1) || = iV], (9) 
1 - P N )a(t - l)B t Q*\l< 

hence, 



f t = (l - p„)a(r - l)«iD-«|« < N2 - (r - D^C^" 1 ]; 



|f + (r - i) || = a(r - 1)0 - o(r - wp*. uo) 

The infinite trunk group does not require the analysis of (6) nor the 
boundary conditions (£ = N) of (7), hence (9) becomes 

|2 - (f - i)g -a - D^csr*- 1 ]. (ii) 

Thus, the corresponding equation for the infinite trunk group is 

|?+(r-i> ff = ««•-«»• (12) 

The binomial moments /3,(J, iV) are defined by 

w,m = ±Pj(t,N)( J s y d3) 
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in which (i) designates the binomial function defined by 



(0 = - 



-D---(i-« + D s>1 

(14) 



The binomial moment generating function l(t, w, N) = £,"=0 /3,(J, N)W 
is given by 

J(*, w, N) = g(t, 1 + w, N). (15) 

Thus, the differential equation satisfied by I is 

d l + w ^L = awl- aw(l + w)*Pw. (16) 

The corresponding equation for the infinite trunk group is 

* + .£-«* (17) 

Equation (16) is a linear partial differential equation that can be solved 
by the following device (method of characteristics). Let be a new, 
independent variable and set 

1 = Ke), 

w = w(d), (18) 

t = t(6). 

Then, comparison of 

dl_ . _ dl dt_ , Bl dw ,-gv 

dd~ dtdd + dwdd K ' 

with (16) yields the set of equations 



-jz = aw — aw{\ + w) N P N , 
do 

dw 

dJ = w > 

<fc _ 
dO ' 



(20) 



whose solution for I is then obtained. To exhibit the solution conveni- 
ently, let 

«(*, r) = e- t f t e'a(s)ds ) (21) 
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and 

A = A(0 = a(t, 0) = e- l *a(t), (22) 

where * designates convolution product over (0, t). The solution of 
(16) takes the form 

l(t, w, N) = 1(0, we- 1 , N)e*» 

- w f e™+'-'[l + we r -q N a(T)P N (T } N)dr. (23) 
Jo 

Similarly, the solution of (17) is 

l(t, w, oo ) = 1(0, we- 1 , oo) e A«- > (24) 

To obtain the binomial moments themselves, it is convenient to intro- 
duce the Volterra operator K„ 

K.f = ['K.it, r)/(r)dr, (25) 

Jo 

denned by the kernel 

K.(t, r) - if £t N ) *-(-i><«-')a(T). (26) 

J =0 J '■ \ i> 1 J / 

Since the Laguerre polynomial L£ a) ( — x) 10 is given by 

the kernel K,(t, t) may be written more compactly; thus, 

K.(t, r) = a(r)e-*<<-'>LJV +1) (-«e'- T ). (28) 

The following theorems may now be stated. 

Theorem 1: The binomial moments, @„(t, <x>), for the infinite trunk group 
are given by 

m, «o = i /wo, •)#-<-*■ £■ 

Proof: The coefficient of w' in the expansion of the right-hand side 
of (24) yields the result. 

Theorem 2: The binomial moments, &,(t, N), for the finite trunk group 
satisfy 

P.(t,N) =(J.(t, oo) - K,P N . 

Proof: The coefficient of W in the expansion of the right-hand side of 
(23) provides the required result. 
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Corollary 1: The probability, Pjfit, N), that all trunks are busy satisfies 
the integral equation 

P N (t,N) = M«, °°) ~ KnPn, 
in which 

fcrfc ») = E Pw-tiO, N)e-W-fi'^ 

i =o J I 

Kjv(«, t) = aWe-vv-^L^ii-ae 1 -'). 

Proof: For the finite trunk group (13) shows that 

M«, N) = P N (t, N). 

Hence, the integral equation follows from Theorem 2. The explicit 
expressions for 0n(1, °°) and K N {t, r) are obtained from Theorem 1 and 
(28), respectively. 

The special case of all trunks free initially leads to a somewhat 
simpler integral equation for Pn(t, N). This is given in the following 
corollary. 

Corollary 2: When all trunks are initially free, P^it, N) satisfies 
PAt, N) = ^ - K N P N . 

Proof: The initial probability distribution P;(0, N) has the form 
P,(0, N) = 1 j = 0, 

= j > o. (29) 

Hence, the binomial moments satisfy 

0.(0, tf) = 1 « = 0, 

= « > o. (30) 

The equation for /3jv(^, °°) given in the first corollary now yields 

M«, •> - jn' (31) 

The result of the corollary follows. 

The probability Pat (2, N) corresponding to all trunks initially busy 
is called the recovery function; it satisfies the following integral 
equation. 

Corollary 3 : When all trunks are initially busy, one has 
P N (t, N) = er"*L N (-e'A) - K N P N . 
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Proof: We have 

Pj(0, N) = ^ j < N, 

= 1 j = N. (32) 

Hence, the required binomial moments are 

0.(O,N) = ( N 8 )- (33) 

The result follows from the equation for 0n(t, «>) and from (27). 

A noteworthy case occurs when a is constant, then the integral 
equation for Pn(1, N) becomes of convolution type. Thus, let 

K N (t) = ae-x'Ltil^-aie* - 1)]. (34) 

Then one has Corollary 4. 

Corollary 4- When the offered load a is constant, PnU, N) satisfies 
P*(t,N) =0 N (t, oo) - Kh*Pn. 

Proof: Use of (21) and Ks{t, r) as given in Corollary 1 yields the result. 
For constant a, an equilibrium distribution P,-(«, N) exists. 2 Let 



ftrW-Sfr (35) 



Then 



^"•^-JTEW O^^iV. (36) 

The notation B(N, a) is used for the blocking probability Pn( x , N) 
and is referred to as Erlang's loss formula. 6 Corresponding to the 
equilibrium distribution, one has the binomial moments /3,(0, N) and, 
hence, the moments for the infinite trunk group given in Theorem 1. 
These moments will be denoted by fi e s (t, a). It may be noted that 

lim/3J(«, a) = ~- (37) 

The integral of K N (t) that is useful in error analyses may now be 
easily obtained. 

Theorem 8: When the offered load is constant, we have 

h KN{T)dT B{N~d) ~ l ' 

f K N { T )dr = S N (a) - 1, 
Jo 

in which K n (t) is given in (34). 
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Proof: Since B(N, a) is the equilibrium solution of the integral equa- 
tion of Corollary 4, and since fir(t, a) corresponds to the equilibrium 
state, the solution of the integral equation is constant and equal to 
B(N, a); thus, 

B(N, a) + K^BiN, a) - &r(t, a). 

This immediately implies the first equation of the theorem. The second 
follows on considering t — ><», and using (36) and (37). 

It may be useful to observe that the positive character of the general 
operator K N in Corollary 1 immediately supplies the inequalities 

M«, •) - K N p N (t, oo ) < P N (t, N) < M«, °°)- (38) 

A Neumann-series solution of the integral equation of Corollary 1 is 
discussed in Appendix A. Higher-order inequalities of type (38) are 
also given. 

III. THE INFINITE TRUNK SOLUTION 

It will be convenient to express the solution for P x (t, <») in terms of 
Poisson-Charlier polynomials 11 whose relevant properties are dis- 
cussed in Appendix B. The probability distribution of the number of 
busy trunks for the infinite trunk case was considered by Khintchine 2 
and, for constant offered load, by Karlin and McGregor. 12 Theorem 4 
presents a rapidly convergent form of the solution in terms of Poisson- 
Charlier polynomials valid for any initial state. This solution will 
be the main tool for the construction of approximations to distribu- 
tions in the finite trunk case. 

From (15), let 1(0, w, ») be given by 

1(0,10. oo) = L/3y(0, 00)1*/, (39) 

then the binomial moment generating function for the infinite trunk 
case is, from (24), 

l(t, u>, oo) = e A - E fr(0, oo)r*w'. (40) 

3=0 

The mean, /x, of this distribution is the coefficient of w ; hence, 

n = A + /xoe-S (41) 

in which no is the mean of the initial distribution. One may now state 
Theorem 4. 
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Theorem 4- The probability distribution, P x (t, «), of the number of busy 
trunks in an infinite trunk group has the convergent representation 

P x (t, «) = *(*, m) t r*d&fa m), 

j=0 



~"ii v ! 
provided that 



d,=- £ ^rf 1 M'/S,w(0, «), 



1(0, • - 1, oo) - £ fr(0, »)( z - l)i = £ P,(0, oc) 2 i 

converges for \z\ < r(r > 2). 

Proof: Use of (40) in (101) which the choice X = n as given in (41) 
yields 



Hence, 

and, from (95), 

Let 



C(w) = e-w>»« ' £ fr(0, »)«-»*«?'. (42) 

1-0 



cy = r* £ ^f-' Mft-»(O f • ) (43) 



P.ft «) - Hx, m) £ <*GM*, M). (44) 

y=o 



then 



rfy = L o ^p Mo%-,(0, co), (45) 

Cj = e-ttdj, (46) 

and the formula of the theorem follows. A theorem of Uspensky 13 
states that the general representation of (95) is valid in the sense that 
the series converges to F(x) if the series £"= F(x)z x has radius of 
convergence greater than 2. Since 

l(t, 2-1, ») - e- A+A *Z[0, (z - l)e- ( , * ], (47) 

the radius of convergence of l(t, z — 1, co) is greater than 2 by the 
condition on 1(0, z — 1, oo) stated in the theorem. Hence, by 
Uspensky's theorem, the representation of (44) is valid for all t ^ 0. 
A truncation error estimate for the series representation of Theorem 
4 is given in Theorem 5. 
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Theorem 5 : 

\P,(t, °°) - +(*, m) - *(a, m) E e-HjQjix, /*) 

i=2 



(0\A+1 g-(A+l)t 

«) 1 - (2/^ ^°^ "^ 



fc £ 1. fl > u> >0, * > In - , 

in which R is the radius of convergence of 1(0, w, <»). 

Proof: Since 

\f(x, /*)<?*(*, m) I = l* {/) (s, m) I ^ 2', j ^ 0, (48) 

we have 

£ r-*|4i| |*(*, M)Gi(*, m)I ^ E e-^|d,-|. (49) 

Also, from (45), 

< V ? 

,1 



dy| ^ E -,/3i-»(0, ») ^ e»°"l(0, w, oo)ur'. (50) 



j=0 v 

Thus, 

E e- y '|dy| Itffa n)Gj(x, ft) \ ^ e»°H(0, w, oo) £ <r*2%-' (51) 

j>k i>* 

5) 1 - (2/«)«- ' (52) 

The conditions of the theorem ensure the convergence of Z(0, w, ») 
and of the series of (51). 

The corollaries below follow directly under the conditions of Theorem 
5. 



Corollary 1 : 



P.(«, «) = *(*, m) E rHflfa n) + 0(«r<*+"9. 

j=o 

Corollary 2: 

\P,(t, .) - *(*, ,)| g | j _ g^p «~t(0, », «). 

Corollary 8: 

PS, •) = f(*. M) + 0(e-«). 

Thus, the distribution quickly becomes nearly Poisson with the time 
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variable parameter n, regardless of the initial distribution. In fact, if the 
initial distribution is Poisson, one might anticipate P x (t, °°) to remain 
Poisson for all/ ^ 0. This is asserted by Theorem 6. 

Theorem 6: If P x (0, °o) is Poisson with parameter n Q , then 

Proof: The binomial moments, /3„(0, <x>), are 

0.(0, °°)=7;, s^O. (53) 

It follows, from (45), that 

d, = ff£ o (f)(-l)" = 0, i>0. (54) 

The result is now obtained from Theorem 3. 

IV. APPROXIMATIONS 

The Neumann series (75) is an explicit solution of the integral 
equation of blocking given in Corollary 1, Theorem 2. For constant 
offered load a, the resolvent kernel solution (108) and Theorem 11 are 
also available; however, especially when N is large, these solutions 
are not convenient. It is therefore important to have approximations 
that lend themselves to computation for large N in a sufficiently 
convenient manner. Three such approximations have been developed, 
namely: the "modified offered load" approximation that is useful and 
fairly accurate in the general case, that is, for time variable offered 
load, the "scaling" approximation, and the "Laplace transform" 
approximation, which are applicable only to transient phenomena 
under constant offered load. The scaling approximation is also fairly 
accurate and does not require factorization of polynomials. The Laplace 
transform approximation consists in fact of an infinite set of approxi- 
mations of arbitrarily high accuracy. It usually requires finding a 
single root — the so-called dominant root — of an appropriate poly- 
nomial. The scaling and Laplace transform approximations are dis- 
cussed in Appendix D. Appendix C provides approximations for the 
required dominant root. The modified offered load approximation 
is presented below. 

Let Pi, x (t, N) be the probability that the iV-trunk group started 
from state i at time and proceeded to state x at time t, then Pi, x (t, a> ) 
may be computed from Theorem 4 using 

&(0, °°) = Q)- (55) 
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An approximation, P ilX (t, N), for P ilX {t, N) is given by 

PUt,N)= /'-<'■*> • (56) 

E Pi At, «) 

..=0 

This approximation is suggested by the following considerations. One 
has 

lim Pi, x (t, N) = Pi At, »)• (57) 

Hence, the approximation should be accurate even for time-varying 
offered load when N is large. Furthermore, when a is constant, since 
lim*-,,*, n = a, one has 

Hm PUt, N) = -nf-TZ-, , (58) 

which, as previously indicated, is the exact equilibrium distribution; 
hence, the approximation should be accurate when t is large even for 
time- varying a(t), provided d(t) is small. Since, by the law of total 
probability, 

P x (t, N) = £ P,(0, N)P<At, N), (59) 

i=0 

one can construct an approximation, P x (t, N), to P x (t, N) by use of 
Pi.S, N); thus, 

P x (t, N) = £ P,(0, N)PUt, N). (60) 

To facilitate the use of (60), Theorem 7 expresses P ilX (t, ■») in finite 

form. 



Theorem 7: 

x\ 



PUt, *) = (!- 0*e- A ^^- (?,[*, -(«' - DA]. 



Proof: From Theorem 4 and (55), we have 

Pi At, «) - *(*, m) £ er'-'djGiix, M ), (61) 

*-£*=£«•(,' )• (62) 

Comparison of (62) with (90) shows that 



<i/ = Tl <?*(*,*). (63) 
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Let 



then 



Cy-^p-'flW^O, (64) 

C(w)= ±^i Gj {i,i)w*, (65) 

which, by comparison with (89), may be rewritten 

C(w) = e-'-'^l + e-'w)\ (66) 

Let 

B(w) = e* w C(w) = e Kw {\ + erhoY (67) 

and 

p(«) - B(z - 1) - e- A e A «(l - er* + rt)«. (68) 

Then 

Pi Ah oo) 

= (1 - r-0«(e* - D-e- A £ ( i ) [(e< ~, 1)A:|y - (69) 

„=o \ £ — v / vl 

Hence, comparison of (69) with (90) finally yields 

PiAt, <*>) = (1 - <r')vr A ^j^ (?,[*, -(e' - 1)A]. (70) 

Immediate corollaries are the following : 
Corollary 1 : 

P.(*, co) = e -A i=^! £ Pi(0> oo) (1 - ryojfr -(e* - 1)A]. 

X, \ i =o 

Corollary 2: 

P ix(t) N) = Mh ~( et - X M 

E o (x \/v !)(-A)« (?,[*, - (e« - 1)A] 

Of particular interest are the functions Po.xit, N) and Px.Nit, N) ; 
the first describes the progression of the system from initially empty to 
blocked, and the second describes the recovery of the system from an 
initially blocked condition to the blocked condition again. The latter 
function is called the "recovery function." 7 The following formulas 
are obtained from Corollary 2. 

ftjrft N) = B(N, M ). (71) 

The general principle of approximation employed, namely, the re- 
normalization of an appropriate solution for the infinite trunk group 
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case, allows one to state, by use of Theorem 6, that, whenever one starts 
from an approximately Poisson initial state, an approximation to 
Ps(t, N) is 

PnU, N) = B(N, M ). (72) 

This approximation was known to C. Palm. 1 The parameter n is re- 
garded as a modified offered load. 

The recovery function approximation obtained from Corollary 2 is 

GnIN, - (e' - 1)A] 



Pw, N (t, N) = -$ — 

£ (#!/,!)(- A) '-"(?,[#, -(e' - 1)A] 



(73) 



■ =o 



V. NUMERICAL EXAMPLES 

For the purpose of providing some idea of the accuracy of the ap- 
proximations developed in Section IV and Appendix D, curves were 
drawn up comparing exact and approximate solutions for a group of 
ten trunks. These curves illustrate nonstationary behavior. Figure 1 
shows the scaling and modified offered load approximations for a step 
input problem in which a = 7 erlangs is offered to an initially empty 
group. The scaling approximation of (166) was used, and (71) was used 
for the modified offered load approximation. Apparently, for this situa- 
tion, the scaling approximation is somewhat more accurate. 
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Fig. 2 — Recovery function — Laplace transform and modified offered load approx- 
imations. 

The recovery function approximations of (73) and (187) are com- 
pared to the exact solution of (185) in Fig. 2. The approximations are 
correct at the extremes t = 0, t = °o , and track the exact curve reason- 
ably well. The approximation of (73) is more accurate initially and is, 
of course, also applicable when the offered load is time- variable ; how- 
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Fig. 3 — Pulse response — modified offered load approximation. 
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Fig. 4 — All trunks empty initially — modified offered load approximation. 
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Fig. 5 — All trunks empty initially — modified offered load approximation. 
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Fig. 7 — All trunks empty initially — modified offered load approximation. 
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ever, the approximation of (187) is simply the case n = of Theorem 
16. Considerable enhancement of accuracy may, for example, be ob- 
tained by using n = 1 or even higher values of n. 

The modified offered load approximation, (72), is compared to a 
pulse response in Fig. 3. It is seen that, despite the rapid variation of 
a(t) (as high as 14 erlangs/call duration) and the large range of prob- 
ability values, the approximation well imitates the course of the 

response. 

Figures 4, 5, 6, and 7 illustrate the modified offered load approxi- 
mation applied to sinusoidal inputs. In all cases, the trunk group was 
initially empty. Figure 4 shows the response, starting from t = 0, to 
a = 7 + 3 sin 4tt«. This may be considered to have wide excursions 
compared to the constant term 7, and rapid oscillations, i.e., the period, 
T, is 5. The exact curve is seen to be well imitated by the 
approximation. 

One may consider the total error to consist of two components, an 
evanescent part arising from the specific initial state and a component 
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Fig. 8 — Upper bound approximation to dominant root. 
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resulting from the rate of change of the offered load itself, that is, a 
function of d(t) which vanishes when d(t) = 0. Figures 5, 6, and 7 
are intended to illustrate the latter component above ; hence, the time 
scale starts at five. The periods of 40, 20, and 4, respectively, were 
chosen to reflect the effect of d(t) on the approximations. To provide 
clearer comparison, the probability scales of the graphs have been 
expanded. 

The lower bound of Theorem 13 provides an upper bound on the 
dominant root. A comparison with the exact values for an eight-trunk 
group taken from Benes 7 is given in Fig. 8. 

VI. NEEDED INVESTIGATIONS 

Much work remains to be done to provide a satisfying and viable tool 
for fully available trunk group analyses. One may mention error estima- 
tion of the approximations suggested in Section IV and Appendix D, 
the investigation of new approximations, such as studying the con- 
sequences of using a refined scaling approximation or an improved 
modified offered load, and the study of (/ + Kn)- 1 in the Banach 
space X, for variable a(t). This, in turn, would permit new approxi- 
mations to be constructed and would provide improved means of error 
investigation. Relatively little is known about the behavior of the 
zeros of Poisson-Charlier polynomials, especially in the present 
context, as functions of a, N. It is hoped this paper will provide an 
impetus for further investigation. 
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APPENDIX A 
Neumann-Series Solution 

The Neumann-series solution of the integral equation for blocking, 

P N (t, N) = fait, co) - K N P N , (74) 

is 

P N (t, N) = N (t, co) - K N p N (t, oo) + K 2 N N (t, co) - . . ., (75) 

which, of course, is convergent for all t. The positivity of K N implies 
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the system of inequalities 

E (-iyK } N p N (t, °o) < p N (t, N) < £ {-iyK> N p N {t, oo), k ^ o, 

(76) 

which generalize (38). 

Truncation error estimates for (75) will now be obtained. The in- 
equalities (76) yield 



-K 2 N k+1 p N (t, oo) < P N (t, N)-Y, (-D'KtfNd, co) < 0, (77) 



z 



2fc+l 



< P N (t, N) - £ (-1)«&M<, oo ) < K% k+1 N (t, co). (78) 

i-o 

Hence, it is only necessary to bound X^+^wfo oo). Let 

o(0 ^ o, (79) 

then, from (27) and (34), 

K N (t, r;a) ^ K N (t - r;a) < Ae-«-*\ (80) 

A = aL&Ui-a). (81) 

The dependence on a is explicitly shown in (80) . One similarly obtains 

fcrft - 5 a) ^ jS = £ /WO, AT) f-! , tfe 0. (82) 

One may now state Theorem 8. 
Theorem 8: 

r\ 
Proof: One has 



&(*;«) ^ Ae-'. (83) 



Hence, 

K NtT (t; a) ^ A r (r _ 1) , 



(84) 



in which K N , r {t; &) is the r-fold convolution of K N (t; a) with itself. 
Convolving this with fi N (t, ») finally yields 



K N Af,*)*Mt, »)^^e-'. (85) 
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APPENDIX B 

Poisson-Charlier Polynomials 

Some properties of Poisson-Charlier polynomials 11 are developed in 
this appendix, especially with a view to convenient representation of 
functions by series expansion. Let 

fix, X) - e~ x ^ , s-0,1,2,--.. (86) 

Then the polynomials Gj(x, X) are denned by 

£?#(*, X) = V j) (x, X) = Gj(x, \)f(x, X). (87) 

The Taylor expansion 

^(x,X + D = jC o j[*»(x, X) (88) 

yields the generating function 

e ''( 1 + i)' = io°' (x - x)t f\- m 

Thus, explicit formulas for Gj(x, X) are 

G,(*,x)=f; £(-«•(.• )^ 

X J „=o Vj — j» / v! 

-i^^O'^O)- (90) 

The first few polynomials are 



Gofa X) = 1, 

:* - x), 

(91) 



Gi(x, X) = i (« - X), 



G 2 (a:, X) = i [a; 2 - (2X + l)x + X 2 ], 
£,(&, X) = i [> 3 - 3(X + l)x 2 + (3X 2 + 3X + 2)x - X 3 ]. 
A recurrence relation derived from (89) is 

G j+1 (x, X) = x - * - X G,(x, \) - 1 Qj-iix, X). (92) 
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The following inner product is denned for functions f(x), g(x) of the 
discrete variable x : 

(/, 9)=t H*, »f(x)g(x). (93) 

i=0 

The Poisson-Charlier polynomials are orthogonal with respect to this 
inner product 11 

(G h G k ) = 0, j*k 

(G h G s ) = i}- 

Accordingly, the coefficients Cj in the expansion of a function f(x) 
in the form 

f(x) = f(x, X) £ ofifa X) (95) 



are given by 

/» . = — 



c t = i! £ Gsix, \)fix). (96) 

j! x=o 



For the purpose of the present investigation, a more convenient mode 
of determining cy is achieved by obtaining their generating function; 
that is, 



Ciw) = £ c/w*. (97) 

j=0 



From (96), one has 



C(w)= £ f(x) £<**f' <?,(*, A). (98) 

i =0 j -o J I 



Hence, from (89), 



C(«0 = «-*" £ /(*)(1 + »)'. (99) 

i=0 



Let 

Biw) = £ fix)il + w) x . (100) 



i=0 



Then Biw) is the binomial moment generating function of fix) and 

dm) = e-^Biw). (101) 

From (101), the first few coefficients Cy are obtained in terms of the 
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corresponding binomial moments ft of f(x) ; thus, 

Co = 0o, 

d = 0i - X/3 , 

c 2 = ft - X/3, + §X 2 /3 , ( 2) 

c 3 = /3 3 - X/3 2 + £A 8 0i - |X 3 /3 . 

A useful choice of the parameter X is suggested by (102), namely, 

X = /V/3o, (103) 

which yields ci = and 

*-*-£■ 



C3 = ^3 - ^ + ^ 



0^2 , ft 

0o 2/3g 



(104) 



For a probability distribution, one has /3 = 1, and the choice (103) 
for X implies that X is equated to the mean of the distribution. In this 
case, one has 

C2 = K«r 2 -M), 

C 3 = *(« - 3<r* + 2 M ), ( U5} 

in which n is the mean, a 2 the variance, and a the third moment about 
the mean of the distribution. 

APPENDIX C 

Constant Offered Load — Dominant Roots 

In this appendix, the integral equation for the constant a case, 
namely, 

P N (t, N) = M«, °°) ~ K N *P N (t, N), 

K N (t) = ae- Ar 'L,y_ 1 [-a(e' - 1)], (106) 

a N (t, oo) . e->n £ ^_ y (0, AT) [a(e< ~ 1)] ' , 
7 =o .? ! 

is studied. In fact, the somewhat more general equation, 

f(t) + K N *f(t) = g(t), (107) 

is resolved. This presents a considerable advantage over the solution 
of (106), since the errors of approximations to P N {t, N) satisfy (107), 
and hence may be studied by means of Theorem 10 and its corollary. 
Solutions for blocking have been obtained in the literature, 78 but do 
not provide a means for error analysis. For the practical utilization 
of the solutions, bounds for the exponents occurring in the explicit 
representation of the resolvent kernel will also be obtained. 
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Integral equation theory 14 asserts that a resolvent kernel, Qir(t), 
exists with the property 

/(<) = 0(t) + QN*9(t). (108) 

Laplace transformation will be used to study (107) and (108). 

Theorem 9: The Laplace transforms, K N (s), Qn(s), of K N (t), Qs(i), 
respectively, are 

~ M _ (-l)VG N (-s-l,a) _ 

n , v (« + !)•••(« + #) - (-l) N a N G N (-s- I, a) 
yjrW ~ (-1) VGat(-s - 1, a) 

Proof: One has, from (27), 

w- 1 / N \ X> 

MW-*)- 5 (, + ,)£. < 109 > 

and hence, 

Thus, 

K N (s) = a N £ ( . J , ) ?-! f arC»+-*«(l - e -')'d*. (Ill) 

Letting X = e~ l , one has 

*" (s) - • % ( i + 1 ) $. /.' x * ,+, ~'" 1(1 - x)ldx - (112) 

The integral in (112) is the beta function, B(N + s - j, j + 1). 
Hence, 

^ (s) = a J ?oij + iJ a r(Ar + s + i)- (113) 

One has the following transformations : 

*^-&(i W-<y+-j+i) ' (114) 

^"° (. +1 )- 1 -(, + iV) ^(?)"(«+i)-(H-y-j), (us, 

**»- (h-d-^jo |,( i ) a "' (s+1) - (s+i) ' (116) 

^-w^Hwi^cji^-r)- <117) 
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Thus, the required formula for Kn(s) is obtained from (117) by 
comparison with (90). From (107), one has 





/ + K N *f = g. 


(118) 


Hence, 


i I K N „. 


(119) 




J l+K N s 1+K N *" 


thus, 


CbrO- ,H S \,- 


(120) 



This, together with the formula for K N (s), yields the required expres- 
sion for Qn(s). 

Corollary : 

p,m>, tn - (-^- ^ff-Ti! a) - 

Proof: Use of Theorem 2, Corollary 4, Theorem 9, and eq. (178) with 

fii(0, N) = (j). 

An expression equivalent to this corollary was given by Takacs. 16 
It will be useful now to introduce the Banach space, X, of functions 

f(t) that are bounded and measurable over (0, °o) and normed by 

11/11 -sup |/(0 1 • (121) 

One may now state Theorem 10. 

Theorem 10: The operators I + K N and (I + Kn)- 1 are bounded; 
further, 

||7 + K N \\ = Sx(a), 

mi + k N )-*w = i+r \Q»(t)\dt, 

Jo 

in which the operator norms are those induced by (121). 

Proof: The quantity ||/ + Kn\\ is obtained directly from Theorem 3 
and the formula 

||7 + £jr|| = 1 + J" K„{t)dt. (122) 

Since the polynomials G N (x, a) are orthogonal over (0, <»), it follows 
that the zeros, Pi, • • • , Pn are distinct, real, and positive ; hence, the 
zeros, ri, • • •, tat, of Gn{ — s — 1, a) are distinct, real, and less than 
minus one. The Paley-Wiener theorem 9 applied to Qn(s) now asserts 
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that / + Kn has a bounded inverse with norm C given by 

C = \\(I + K N )-^\\ = 1 +/ o °° \Qx(l)\dt. (123) 

An immediate corollary is the following. 

Corollary: The integral equation 

f + K N *f = g, gGX, 
possesses a solution / €■ X satisfying 

11/11 ^ cm, 

lmi |/(0| ^ CBm |flf(<)|. 

f-»oo t-»oo 

7n particular, if lim t ->*> g(i) = 0, 2/ie?i limt-a, f(t) = 0. 
Proof: The result follows from 

/ = flf + Qn*Q (124) 

and the definition of C. 

The following theorem provides a representation of Qir(t) and an 
estimate of C. 

Theorem 1 1 : 

N ft fa+") at II \r$+v\ 

Qn(1) = E # •"•. C ^ 1 + L ~ ^ 

J=1 n (rt-r,) 3=1 \rj\ J[ \r, - r,\ 

Proof: One has 

(-I)VGat(-s- 1, a) = (s-n)---(s-r^) (125) 

and, hence, from Theorem (9), 

&,(.) . ('+»■••(; + g> - (. - r,)- • •(« - -») . (126) 

* w (a — ri) • • • (s — r*) 

The partial fraction expression for @tf(s) is 

A' 



* n fa + r) j 



Q«(s) = E y firy: (127) 
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Thus, 

N n ixi + v) 

Qn(0 = L V 1 erit - ( 128 ) 

;=1 n (fi-r,) 

v = l 

Also, one has 

at n ki+vi 

[Qi^rCOl ^ L # e" 1 , (129) 

* =1 II |r,--r,| 



and, accordingly, 



A' 



,„ * nir,-+H 

C = l + / IQ^WIA^l+E— ^ (130) 

; ° J=1 \rA n |ry-r,| 

To make the results of Theorem 1 1 accessible for estimation, partic- 
ularly in numerical applications, and for the Laplace transform ap- 
proximation of Appendix D, bounds for the roots, ry, will be obtained. 
The generating function, g(z), for the equilibrium probability distribu- 
tion, Pj(°°, N), of (36) may be written as 

Thus, the mean, m, and variance, a 2 , are 

m = a[l - B(N, a)], (132) 

ff t = m - (N - m )aB(N, a). (133) 

We now have Theorem 12. 

Theorem 12: (Benei) ri > — (m/<r 2 ). 

To obtain further bounds, the following lemmas are needed. 
Lemma 1: pn ^ ■ • • ^ pi > 0, 

Pi + • • • + PN = Si, Pi + * ' • + P% = S 2 , Pi' ■ -PN = D, 

then 

< . . . < 



^ 8l + V(JV- l)(ATs 2 -s 2 ) 



2?i which p is the small positive root of 

P ( Sl - p)»-* = D(N - !)*-». 
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These bounds are sharp. 

Proof: The equations 

(134) 

(135) 

(136) 
(137) 

(138) 

over all allowable sequences. Then the sum (137) is minimum when p# 
is replaced by p* N . This occurs, however, only when 

pi = ■ • • = pjv-i = ft _ i " (139) 

Thus, from (137), 

(tf-l)(^3T) 8 = s *-pft (140) 

The solution of (140) for p* N is 





Pi + • • • + PN = Si, 




Pi + • • • + Pn = s 2 


may be written 






Pi + • ■ ■ + PJV-I = Si — Pat, 




Pi + • • • + Pat-i = S 2 — p%- 


Let 






p* N = max pn 



, Si + V(AT-l)(^ S2 - S f) 
Pat = jj (1*1; 

This proves the upper bound of the lemma. The inequality is attained 
for the vector (pi, • • •, par) denned by (139) and (141). 
Similarly, one may write 

P2 + • ■ • + pat = s x — pi, (142) 

P2---PAT = D/pi. (143) 

Let 

p? = minpi (144) 

over all allowable sequences. Then the product (143) is maximum when 
Pi is replaced by p*. This occurs only when 

* 

p 2 = • • • = PN = ATZTf ' ( 145 ) 

Thus, 

[ w^t }"" = D/ "' (146) 

and the equation of the lemma defining p has been established. The 
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inequality is attained for the vector (pi, • • •, />#) denned by (145) and 
(146). 



Lemma 2: />\ ^ 



^ Pi > 0, 

1 



s_i = - + 
Pi 



+ 



D = 



PN 



Pi' * 'PN- 



> 



> 



> 



Then 

PN ^ • • • ^ P 1 ^ P, 

in which p is the small positive root of 



dL-i-^Y" 1 = P (N - i)»-K 



The bound is sharp. 






Proof: One may write 






pi 


. 1 1 
• • H = s_i , 

PAT Pi 


(147) 




p 2 . . . pAr = jD/pi. 


(148) 


Let 


pt = min pi 


(149) 



over all allowable sequences. Then the product (148) is maximum 
when pi is replaced by p*. This occurs when 



N - 1 



Pi = • • • — PN — 



S-l — 



Pi 



Thus, 



r n - 


- 1 ] 


S-i - 


1 

* 

Pi. 



= D/ P \ 



(150) 



(151) 



and the equation of the lemma is established. The inequality is at- 
tained for the vector (pi, • • •, pat) defined by (150) and (151). 

The application of the lemmas to the polynomials Gn(x, a) is ac- 
complished by identifying pi, ■ ■ • , pn with its zeros. For this purpose, 
the form of Gn(x, a) given in (90) will be recast, by the help of the 
Stirling numbeis of first kind, 11 into standard form; that is, 



G N (x, a) = Y. a N - m x m . 

m=0 



(152) 
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The Stirling numbers, Sf, are denned by 

ff (x - v) = t S?* m - (153) 

►-0 m=l 

Thus, one has 

G„(x,a) = (-1)" + £ x» N £ g (-l)'a'-»( N )sS-.. (154) 

The sums s x , s 2 are accordingly given by 

s 1 = (^ + aN, (155) 

s 2 = sf - e(^) - (6a + 4) (£) - 2a 2 (^)- (156) 

The reciprocal polynomial, that is, the polynomial whose zeros are 
pr 1 , • • •, Pn 1 is given by 

X N G N (x~\ a)G N (x, a) = £ a M a;"\ (157) 

m— 

Hence, the analogous quantities s_i and s_ 2 = p^ 2 + ■ ■ ■ + pn 2 are 
given by 

s_i = £ - JV<»a->, (158) 

s_ 2 = s 2 _, - 2 £ i N^ar" J 1 -. , (159) 

in which 

AT<°> = 1, iV» = iV(iV - 1) • • • (N - v + 1), v > 0. (160) 
The upper bound of Lemma 1 now establishes Theorem 13. 

Theorem 18: 

N 



S_!+ V(tf - 1) (#8-, -«»_!) 



^ Pi < • • • < PAT 



8l + V(iy- i)(JVs 2 - s f) 



Proof: The upper bound is immediate. The lower bound is obtained 
by applying the upper bound of Lemma 1 to the reciprocal equation. 
A numerical illustration of Theorem 13 is provided by the zeros of 
Gio(x, 7) used in obtaining the recovery function plotted in Fig. 2. 
They are 0.332811, 2.05847, 4.06653, 6.31227, 8.81308, 11.5197, 
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14.6407, 18.0255, 22.0872, and 27.1438. The lower and upper bounds 
given by the theorem are 0.32964 and 36.82292, respectively. It appears 
that the lower bound may well be usable as an approximation to pi. 
The accuracy of this when used to approximate —r\= 1 + pi is 
illustrated in Fig. 8. The exact values of —r\ are taken from Bene§. 7 
This provides an upper bound for ri which, together with the available 
lower bounds, is useful for error investigations. 

APPENDIX D 

Approximations — Constant Offered Load 

It was suggested by S. Horing that an appropriate scaling between 
Pjv(£, N) and Pk(t, k) may exist; that is, a function F[_Pff(t, N), 
Pk(t, fc)] may exist, which would be approximately independent of t 
and which may, therefore, permit the approximate determination of 
PnU, N) in terms of Pi(t, 1), for example, thus permitting large trunk 
groups to be studied in terms of the behavior of small ones. Since it is 
feasible to use (108) for small trunk groups, this would constitute an 
approximation of the solution of (106) for large trunk groups. 

Consider the following 

B(kl, a) Wk kl Si(a/k) k nfin 

B(l,a/k) k (kl)\ S k i(a) ' K } 

The ratio Si(a/k) k /S k i(a) is approximately independent of a since 
Si(a/k) k = Ski(a) = e° for fc large. Thus, 

B(kl, a) (162) 



B(l, a/k)> 



is approximately independent of a. It would seem, therefore, that the 
ratio (162) is approximately a time invariant of (106), especially for 
large t ; that is, the function 

P k i(t, kl; a) nfi oN 

Pi(t,l;a/k) k UM; 

is approximately equal to the ratio (162) when t is large; thus, 

«'i«' s W»[w]' (164) 

We have, from (107) and Theorem 11, 

Pl(U;a) = ^Tl u " e " (0+1) ' } ' (165) 

when the trunk is initially empty. Hence, from (164) with I = 1, 
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k = N, 

P N (t, N)^B(N, a){l - e-K-W+i]'}* ( 166) 

for the case when all N trunks are initially empty. 

It is desired now to extend the approximate solution (166) for 
application to any initial condition. Define q^ (t) by 

q N (t) = B(N, a){l - erlMm+i}t}if m (167 ) 

Then a convolution operator with kernel Ln (i) will be constructed so 
that one has 

q N +L N *Q N =j^, (168) 

exactly. The operator Ln will then be taken to approximate the opera- 
tor K N of (106). Accordingly, an approximation P N to Pn correspond- 
ing to initial states other than all trunks empty is defined by the 
equation 

P N + L N *P N = fa. (169) 

Theorem 14: The Laplace transform, Ln(s), of the approximating 
L N (t) is 

/ n \ r ( S ) r ( N + / n /A/Ul + 

l+W£ + lW(a) ^ (a/ / )+ -+ 

Proof: One has 

q N = B(N, a) f°e-*{l - e-K./w+ni}^ (170) 

Jo 

Let u = Z(a/N) + 1J, then 

* N = WN)'+ 1 lo e - [ " (alN+1)]u ( 1 ~ e ~ u ) Ndu - ( 171 ) 

The substitution x = e~ u yields 

fr = (a/W + l / § 1 * l-/<- '* H), - 1 (l " *)*<**; ( 172 > 

thus, 

«* - (oT^TR r r y+ « +1 1 ' (173) 

L + (a/iV) + l + 1 J 
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One also has 

f^-^- rff+i) ' (174) 

The Laplace transform of each term of (168) now yields 

frd+fi. -gffiL^ . (175 ) 

Use of (173) yields the result of the theorem. 

It may be noted that Ln has an impulsive component at t = whose 
value is 

(£ + l)"V(a)-l. (176) 

For large iV, this is nearly zero. The effect of this is to create an 
error at t = 0; that is, P„(0, N) * P N (0, N) if /MO, N) ^ 0. The 
larger N is (for fixed a), the smaller the discrepancy. 
Equation (169) is studied in Theorem 15. 



Theorem 15: The solution of (169) is 



p » = WW+i B(N > a) 



[ (a/N) + 1 J 



r(s)r \n + 



(a/N) + 



I-] 



• £ 0,(0, N)arV(8 + j), 



Proof: From Corollary 1, Theorem 2, one has 



hence, 



= £ /WO, tf) ?. «rc*-j>i(i _ e-«)i ; (177) 



^=|„ ft(0 '^-'f(^¥T-,T < 178 > 



Transforming the terms of (169) yields 

P N (l + Lat) = fc, (179) 

and hence the result of the theorem is obtained from Theorem 14 and 
eq. (178). The inversion of the transform, Pn, by use of differentiation 
follows on use of (173) in Pn. 
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Another method of approximation useful for constant offered load is 
based on an approximate inversion of the Laplace transform. Let 



J( s ) = f M e -™f(u)du, s > 0; (180) 



then 



(-!)• 



n! 



s n+l 



/(n)( s ) = f_^ / e -™u n f(u)du. (181) 

n ! Jo 

The function (s n+1 /n \)e~ tu u n is a probability density function for 
any s > 0, n ^ whose mean is (n + l)/s and variance (n + l)/s 2 . 
Letting s = (n -f- 1)/% the mean and variance are t and t 2 /(n + 1), 
respectively; hence, Korovkin's theorem on sequences of positive 
functionals 16 establishes the Widder inversion formula 9 : 



lim i-4^s n+1 / (n) («) 



- f(t). (182) 

»=(n+l)/< 



The above discussion forms the basis for Theorem 16. 

Theorem 16: For e > 0, let the transform of e lt f(t), namely, f(s - e), 
exist for s > 0, and let e tl f(t) be convex in t > 0, then 



fit) 



( — 1) n 

g e -«t i_LL S n+l/(n)( s _ c ) 



n ^ 0, t> 0. 

«=(n+l)/< 



Proof: Jensen's inequality applied to (181) in the form 

(~V" s n+1 / (n) (s - e ) = s -l- ( K e -° u u n e'*f(u)du (183) 

n! n\ Jo 

establishes the theorem. By virtue of (182), when similarly modified 
for the function e lt f{t), the dexter of the inequality always provides 
an approximation to f{t) even when e tl f(t) is not convex. 

Corollary: f(s — e) exists for s > 0, and e tl f(t) is convex in t > 
implies 

Kt) g i r-/ ()-<). t>0. 

Proof: The case n — of Theorem 16. 

If f(s) should have a dominant pole, it is usually advantageous to 
choose e equal to the negative of that pole. 

The above corollary will now be applied to obtaining an inequality 
for the recovery function. The corollary to Theorem 9 shows that 

sGn{ — s — 1, a) 
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The inversion of Pn,n(s, N) is readily accomplished; the result is 7 

Pn.nH, n) = b(n, a)- e ^ n (i - r^-rV (185) 

It follows from (185) that 

e-^lPN, N (t, N) - B(N, a)] (186) 

is convex for / > and that the corresponding Laplace transform exists 
for s > 0, hence, by the above corollary, one obtains 

Pjw(<, N) < B(N, a) 

+ |_ (1 + ri t)G N l- (1/0 - n - 1, a] 1 + nt J v ; 
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